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An updated formulation of soft diffraction, compatible with s and t channel unitarity, is 
presented. Its consequent general soft scattering features at high energies are explored. 
The critical interplay between theory and data analysis and its consequent implications 
with regards to the theoretical foundations of soft scattering theory are discussed. 



1 Introduction 

The renewed interest in soft scattering and Pomeron physics, which was dormant for many 
years, can be traced to the market demand for reliable estimates of hard diffraction gap sur- 
vival probabilities, notably diffractive Higgs production at the LHC. We note that the present 
vigorous study of this subject is essentially based on sophisticated utilization of relatively old 
theoretical ideas and models, such as Gribov's Reggeon field theory p], Good and Walker (GW) 
decomposition of the proton wave function accounting for low mass diffraction [2] and the eikonal 
approximation |3] which secures the compatibility of the scattering amplitudes with s channel 
unitarity. Compliance with t channel unitarity is associated with multi Pomeron interactions 
(IP enhancement) which are a generalization of Mueller's triple Pomeron mechanism [3], pro- 
vided G^p, the triple Pomeron coupling, is not too small. This mechanism supplements GW 
diffraction with an additional high mass diffraction. 

In this talk I shall discuss the modeling and consequent soft scattering predictions derived 
from the above dynamical considerations. The implied gap survival probabilities will be dis- 
cussed by Gotsman in the following talk. I shall assume a Regge like parametrization in which 
the IP is super critical, i.e. ctp(t) = 1 + Ajp + a'pt, where Ajp > 1. The above IP exchange 
violates s-unitarity at high energies. Recall that implementing s-unitarity is model dependent. 
I shall confine myself to eikonal models which have the virtue of simplicity. 

In the ISR-Tevatron range Otot and a e i are well reproduced by Donnachic-Landshoff (DL) 
iron screened Regge parametrization with Ajp = 0.08 and a' P = 0.25GeV~ 2 . The energy 
dependence of the soft diffractive cross sections (notably a s d) is much milder, implying that 
strong screenings initiated by s-unitarity must be taken into account. As we shall see, the 
interplay between theory and data analysis results in strong constraints on both Ap and a'p 
inputs. This results in profound consequences for the nature of the Pomeron and its QCD 
foundations, providing a unifying interpretation of soft and hard Pomcrons. 
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2 Good- Walker Eikonal Models 



Consider a system of two Fock states, a hadron and a diffractive state ^ d which are or- 
thonormal. The GW mechanism stems from the observation that these states do not diagonalize 
the 2x2 interaction matrix T. Introduce two wave functions and ^ which diagonalize T, 

K£ =< *i *k|T|** **/ >= i4< lfc 4,fe'- (1) 

i.e. the A$ k amplitudes are constructed from the elastic scatterings of vE^ and In this 
representation the observed hadronic states are written 

* /l =a*i+/3*2 *d = +a* 2 , (2) 

where a 2 + /3 2 = 1 . The corresponding unitarity equations are 

ImAf !k (s,b) = \Af tk (s,b)\ 2 + G i i %(s,b), (3) 

where Gfi is the summed probability for all non GW inelastic processes induced by an initial 
(i, fc) state. A general solution of Eq.Q can be written as 

AUs,b)=i(l-e,v(-"^^\Y (4) 

G»( S ,6) = l-exp(-0^( S) &)), (5) 

where ilf k arc arbitrary. In the eikonal approximation flf k are assumed to be real and deter- 
mined by the Born (non screened) input. From Eq.([5|) we deduce that the probability that the 
initial projectiles (i, k) reach the final LRG diffractive interaction unchanged, regardless of the 
initial state re-scatterings, is given by Pf k = exp (— fif k (s, b)). In general, we have to consider 
four possible (i, k) elastic re-scattering options. For initial p-p (or p-p) the two off diagonal 
amplitudes are equal, Af 2 = -i4.fi- The corresponding elastic, SD and DD amplitudes are 

a el (s, b) = i{a 4 ^f,i + 2a 2 l3 2 A s 12 + ^ 4 ^| 2 }, (6) 

a sd (s, b) = ia/3{-a 2 Af A + (a 2 - /3 2 )^f 2 + /3 2 Af 2 }, (7) 
a dd = ia 2 p 2 {Al 1 -2Al 2 + Al 2 }. (8) 

The GW mechanism was originally conceived so as to describe a system of a nucleon plus 
its diffractive N* isobars. Obviously, this simplistic approach is not suitable for high energy 
diffraction where is bounded by 0.05^3, implying a continua of diffractive Fock states. 

Throughout this talk I shall relate to GLMM[S] and KMR[5] models which are conceptually very 
similar, but differ significantly in both their formalism and data analysis. In the present context, 
two procedures were devised to overcome the above difficulty: GLMM lump together all GW 
diffractive states to an effective \D > state, to which we add the non GW IP enhanced high mass 
diffraction. In this approach the GW contribution is very significant and the mass distribution 
is smooth. KMR and LKMR[7] chose to confine GW diffraction to low Mj lff < lOGeV 2 , to 

1 this is an arbitrary bound commonly used. 
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which they add the high mass IP enhanced contribution. In this approximation the bulk of the 
diffractive mass is non GW, and its smoothness at lOGeV 2 is not secured. 

GLMM, KMR and LKMR are multi channel eikonal models in which the initial re-scatterings 
of the incoming projectiles includes also diffractive states. 

Slf lk (s,b) = t% k (8)Tf ik {8,b,...), (9) 

where vff.(s) — <7i5fe(^) Ajp and Tf k are the 6-space profiles. The profiles are external infor- 
mation in as much as, beside their normalization and asymptotic constraints on their behavior, 
they are determined by the data analysis. In GLMM Tf k are given as the ^-transform of a two 
i-poles expression (t = —q 2 ). Setting a'p=0, the profiles are energy independent 
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(1 + q 2 jm 2 ) 2 (1 + q 2 /m 2 k ) 2 
GLMM introduce a small energy dependence 



T s (6; m h m k - a' P = 0) . (10) 



m 2 m 2 {s) = A of,,, \ ■ (11) 

The above parametrization is compatible with the requirements of analyticity/crossing sym- 
metry at large b, pQCD at large q 2 and Regge at small q 2 . For details see Ref.[5]. KMR and 
LKMR use a different parametrization for Tf k which is numerically compatible with GLMM. 
The 3 groups reproduce da e i/dt well in the forward t < 0.5GeV 2 cone. 

Consider a model in which diffraction is exclusively GW. This was recently explicitly consid- 
ered by GLMM and LKMR and non explicitly by KMR. These, as well as earlier GW models, 
fit the (different) elastic sectors of their data bases, obtaining output fitted Ap = 0.10 — 0.12 
with x 2 /d-o.f. < 1.0. The above GW models fail to reproduce the diffractive sectors of their 
data base. This deficiency is traced to the need to add the enhanced 2P high mass contributions. 
This has been done in GLMM and KMR. LKMR enhanced JP input is confined to Eq.(fT] 



3 Multi Pomeron Interactions 

Consider a single diffraction channel p + p — > p + M s( i- Mueller's triple Pomeron mechanism, 
derived from 3 body unitarity, leads to high SD mass which is non GW. In the leading order 

j 3iP i / , \ 2Ap+2a'jpt , M 2 \ AiP 

^^ = ^« o ^(^) b?) • (12) 

The virtue of Eq. (TT2")) is that A p can be determined from either the energy or mass dependences 
of the SD cross sections. This approximation is valid for s >> M 2 d >> m p . 

CDF analysis suggests[8] a relatively large value of G^p- Consequently, we need to consider 
a very large family of multi Pomeron interactions (enhanced JP) which are not included in the 
GW mechanism. As we shall see, this "new" dynamical feature initiates profound differences 
in the calculated values of soft cross sections and non GW diffractive gap survival probabilities 
(soft and hard). Note that this features become significant above the Tevatron energy. 
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GLMM and KMR treatment of IP enhanced interactions stems from Gribov classical Rcggcon 
calculus p] and Kaidalov et al.[9]. Recall that in this context the soft Pomeron is a simple pole 
in the J-planc, while the hard (BFKL) Pomeron is a branch cut. KMR model, which is a par- 
tonic model, derives directly from these sources. Its summation is confined to semi enhanced 
IP diagrams (see Fig. lb). KMR calculations are based on two ad hoc assumptions: 

1) The coupling of a multi IP point vertex nIP — > mP [n + m > 2) is g^ = ^g^nm\ n+m ~ 2 . 
In this notation G$p = Xgiy. Note that in Kaidalov et al. g^ = \gN^ n+m 2 ■ 

2) Most of LHC non GW diffractive reactions of interest are hard. Given a 3IP vertex, G%p is 
unchanged by the interchange of soft and hard Pomerons. This is not self evident. A possible 
support for the above is obtained from GLMM interpretation of the Pomeron (see below). 

As we shall see in the next chapter, the data 
analysis executed by GLMM and KMR converges 
to compatible exceedingly small a'p values and 
high, BFKL like, A P ~ 0.30 - 0.35. The ad- 
justments of these parameters are correlated. In 
non screened Regge model Ap controls the elastic 
cross section energy dependence, while a' F con- 
trols the energy dependence of the (shrinking) 
elastic slope. As a'p gets smaller Ap becomes 
larger initiating stronger screening which compen- 
sates the reduction of a'p, and vice versa. As 
we saw the vanishing value of the fitted a'p was 
explicitly dictated by the b-profiles GLMM and 
KMR chose. Non explicitly this choice initiated 

also the high Ap fitted value. These results have profound implications: 

1) A key observation of GLMM is that the exceedingly small fitted value of a'p implies that 
the "soft" IP is hard enough to be treated perturbatively. Following Gribov we identify the 
correlation between a'p and < p t > , the mean transverse momentum of the partons (actually, 
color dipoles) associated with the IP. < p t >= 1/^Ja'p, from which we deduce that the QCD 
running coupling constant as << 1. Accordingly, we proceed from Gribov's parton model in- 
terpretation to pQCD. GLMM sum over the enhanced diagrams. Technically, we have adopted 
the MPSI procedure[in] in which is reduced to a sequence of triple IP vertexes (Fan dia- 
grams). For details see Ref.[5]. This may pose a problem for the calculation of SD cross section 
for which the lowest order diagram is semi enhanced. To avoid this problem we have added to 
this calculation a term by term summation of the relevant semi enhanced diagrams. 

2) The fitted high value of Ap initiates strong screening which results in a renormalization of 
the Pomeron exchange amplitudes. As a result A C p^ is reduced monotonically with energy. In 




Figure 1: Low order terms of the Pomeron 
Green's function, a) Enhanced, b) Semi- 
enhanced. 



GLMM calculations A^ J (2TeV) ~ 0.070, A e J, r {UTeV) ~ 0.045 and A^ J (60TeV) ~ 0.032. 
KMR results are compatible with ours, see Table 2. The slow decrease of A e J>f raises the 
question if its value may become negative at high enough energies, larger than W = 10 5 GeV 
which is the bound of validity of both GLMM and KMR. In GLMM we have checked that 
A^ / (100TeV") > 0. I am less clear about KMR. The compatibility between GLMM and KMR 
is surprising. As noted, GLMM sum over the enhanced diagrams (Fig. la) while KMR sum 
over the semi enhanced diagrams (Fig. lb). Very intuitively (at the risk of being wrong), it 
seems that GLMM renorrmalize the IP propagator while KMR renorrmalize the IP vertex. A 
complete calculation should, obviously, include both contributions. 
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4 The Interplay Between Theory and Data Analysis 



The goal of the data analysis of interest is to adjust a set of theoretical parameters. To this 
end we construct: 

1) A suitable data base adjusted, with a satisfactory resolution, to fix the parameters. 

2) An adjustment procedure, be it a fit (GLMM, LKMR) or tuning by trial and error (KMR). 
We distinguish between a parameter adjustment which is executed through a reconstruction 
of the entire data base in one step (GLMM), and a sequence of factorizable adjustment steps 
(KMR, LKMR), each adjusting a sub group of the free parameters. 

There is a significant difference between the data analysis carried out by GLMM and KMR. 
This reflects both in the choices of data bases made by the two groups and their adjustment 
procedures. The starting point of both investigations is the observation that a GW model repro- 
duces the elastic data well, but its reproduction of the diffractive sector is deficient. Both groups 
claim to achieve an improved reproduction of their over all data base once the contributions of 
enhanced Pomcron diagrams are included. 

GLMM have constructed a global data base so as to fit the multitude of free parameters. 
It includes atot, cr e i, a s d, &dd and B e i in the ISR-Tevatron range, CDF differential elastic cross 
sections and SD mass distribution were checked for consistency. The conceptual approach of 
KMR and LKMR is completely different. Their data base contains only the measured values 
of duei/dt, which enables to predict (Jtot, and do sd / dtd(M 2 d / s). In my opinion KMR data base 
is too limited to enable a substantiation of their goals. Specifically: 

1) As we saw, the ^-profiles Tf k control the features of da e i/dt which are only weakly coupled to 
the proposed dynamics. Indeed, 6 models, published over the last 10 years (3 KMR+LKMR+2 
GLMM) with different dynamics, i.e., exclusive GW (GLMM, LKMR), GW+zero order IP en- 
hancement (LKMR) and GW+ IP enhancement (GLMM, KMR). The output parameters spread 
over 0.1 < Ajp < 0.55 and < a' P < 0.066. All 6 models reproduce, almost identically, the 
CDF distributions of dcr e i / dt. The unavoidable conclusion is that a reconstruction of da e i / dt on 
it own has no resolution power. The only common ingredient to all 6 models is their compatible 
b-profilcs. These profiles constrain a'p to very small values. This is the key observation leading 
to a pQCD (GLMM) or partonic (KMR) IP interpretation. 

2) The reconstruction of CDF da sd /dtd(M 2 d / s) by GLMM, LKMR and KMR partially sub- 
stantiates the introduction of high mass multi Pomeron interactions. In my opinion, this inves- 
tigation, in its present state state, is unable to provide a decisive verdict on this issue. GLMM 
and KMR sum different sectors of the enhanced diagrams. LKMR take into account only the 
lowest order Mueller's 3IP diagram. Regardless of these differences, the three groups produce 
compatible output results which indicate that CDF data, as is, is not sufficient to differentiate 
between different modes of IP enhanced diagram summations. 

3) To further clarify the experimental limitations, let us recall that CDF define their high mass 

diffraction bound 1 — xl = < 0.15 (the common bound is 0.05). KMR and LKMR define a 
lower bound M s 2 rf > lOGeV 2 which corresponds to 1 — xl = 3 • 10~ 6 . Note that CDF SD mass 
distribution available for analysis corresponds to \ — xl > 2- 10~ 3 . This mass distribution covers 
less than 20% of KMR expected high diffractive mass. GLMM have a completely different clas- 
sification in which the GW low mass diffraction reduces monotonically with no arbitrary upper 
bound and the high mass is defined above lOGeV 2 , identical to KMR. Additional difficulties 
with the analysis of CDF data with 1 — xl > 0.02 is that an arbitrary significant background 
contribution has to be added. It is induced by secondary Regge diagrams such as IPIPR. An 
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added clement of ambiguity is that LKMR with a zero order 3IP calculation is faring as well as 
the high order summations of GLMM and KMR. 





AiP 


a 


a' F 

GeV' 2 


9i 

GeV' 1 


92 

GeV' 1 


TOl 

GeV 


TO2 

GeV 


X 2 /d.o.f. 


GW 


0.120 


0.46 


0.012 


1.27 


3.33 


0.913 


0.98 


0.87 


GW+F - enh. 


0.335 


0.34 


0.010 


5.82 


239.6 


1.54 


3.06 


1.00 



Table 1: Fitted parameters for GLMM GW and GW+JP-enhanced models. 



GLMM fitting procedure aims to reproduce our global data base. A fit with a GW model 
(no JP-enh) provides excellent reproduction of our elastic sector while the the reproduction of 
the diffractive sector is very poor! The repeated fit with a GW+JP-enh model results with 
a very good x 2 - The outputs of both models are presented in Table 1. Checking we note 
that the exceedingly small value of a'p is persistently obtained in both models. The outputs 
of Ajp and 52 change drastically once JP-enh is included. As we shall see this has significant 
consequences for the approach of a e i(s, b) toward the black disc bound. As we have noted, KMR 
and LKMR tune g\ and 52 through a reproduction of da e i / dt and the other parameters through 
a reproduction of da s d/dtd(M 2 d /s). The values they obtain for Ap and a'p are compatible 
with ours. 





Tcvatron 
GLMM KMR(07) KMR(08) 


LHC 

GLMM KMR(07) KMR(08) 


W=10 b GeV 
GLMM KMR(07) KMR(08) 


CTtot 


73.3 74.0 73.7 


92.1 88.0 91.7 


108.0 98.0 108.0 


Pel 


16.3 16.3 16.4 


20.9 20.1 21.5 


24.0 22.9 26.2 


Osd 


9.8 10.9 13.8 


11.8 13.3 19.0 


14.4 15.7 24.2 


Odd 


5.4 7.2 


6.1 13.4 


6.3 17.3 



Table 2: Comparison of GLMM, KMR(07) and KMR(08) cross sections in mb. 



GLMM and KMR high energy Tevatron, LHC and Cosmic Rays predicted cross sections 
are summarized in Table 2. The elastic and total cross section outputs of the two models are 
compatible and, above the Tevatron, significantly lower than those obtained in models with no 
multi-Pomeron contributions. This is a consequence of Ap renormalization due to the enhanced 
IP contributions. GLMM and KMR(07) predicted a s d arc compatible, where KMR(07) arc 
consistently larger by approximately 12%. KMR(08) predicted a s d are considerably larger than 
GLMM and are growing at a faster rate. The difference between KMR(07) and GLMM add 
predictions is even more dramatic, where <Jdd{K M R(07) / OddiGLM M) ~ 3 at W = 10 5 GeV. 
This very large difference is due to KMR large diffractive high mass predictions (see Table 2). 
The recent KMR(08) neglects to mention the high diffractive mass sector of double diffraction 
while showing higher SD cross sections than in KMR(07). 

In my opinion the GLMM and KMR compatible predictions of total and elastic cross sec- 
tions at the LHC and AUGER are of fundamental importance because they arc significantly 
lower than the predicted values based on unitarity models with no IP enhancement. These mea- 
surements may provide a decisive support for the importance of multi-Pomeron interactions at 
high enough energies and, consequently, imply that the growth of the total and elastic cross 
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sections with energy is much more moderate than anticipated. This feature reflects in the slow 
decrease of shown in Sec 3. 



5 The Approach Toward the Black Disc Bound 



The GW base amplitudes of GLMM are Af k , with 6 dependences specified in Eq. (flU|) - Eq. ([TTjl . 
These are the building blocks with which wc construct a e i, a sc i and add (Eq. ([6]) - Eq. ([8])). The 
A? k amplitudes are bounded by the black disc unitarity limit of unity. Checking GLMM fitted 
parameters, presented in Table 1, we observe that <?i and gi, which are comparable in the 
GW model, significantly change in the GW+-/P enhanced model where we obtain g-2 >> gi- 
The implication of of our fitted values of g\ and 172, is that including the diffractive data in 
our global fit forces a large inequality between the three GW Af k components. A§ 2 (s, 6 = 0) 
reaches unity at a very low energy, Af 2 (s,b = 0) reaches unity at approximately W=100GeV 
and Afi(s, 6 = 0) reaches unity at exceedingly high energies, well above LHC. The observation 
that one, or even two, of our Af k (s, 6)=1 does not imply that the elastic scattering amplitude 
has reached the unitarity bound at these (s,6) values. a e i(s,b) reaches the black disc bound 
when, and only when, Af 1 (s, b)=Af 2 (s, b)=A^ ^(s, b)=l, independent of f3. The approach of 
a e i(s,b = 0) toward the black bound depends on the rate of Af 1 (s,b) increase. This rate is 
very slow because of the relative smallness of g\ to which we add that Af 1 increase with energy 
above LHC becomes ever so moderate as a consequence of the renormalization reduction of A p . 
Our results are different from the predictions of most available models, notably KMR, in which 
a e i(s, 6 = 0) reaches unity a few TeV above LHC. Note, though, that GLMM is the only model 
which includes the diffractive along side the elastic data in its data analysis. All models which 
predict saturation of a e /(s, 6 = 0) just above LHC have confined their data analysis exclusively 
to the elastic sector. 

A consequence of the input Clf k being very large at small 6, is that Pf k {s,b) is exceedingly 
small at these small 6 values. As a result, given a diffractive (non screened) input, its output 
(screened) amplitude is peripheral in 6. This is a general feature, common to all eikonal mod- 
els regardless of their b-profiles details. The general behavior indicated above becomes more 
extreme at ultra high energies, where a e i continues to get darker and expand. Consequently, 
the inelastic diffractive channels (soft and hard) becomes more and more peripheral and rela- 
tively smaller when compared with the elastic channel. Given (s,6) at which a e ;(s,6) = 1, the 
corresponding diffractive amplitudes, GW and non GW, vanish. 

The behavior of the ratio Rd = (c e ; + cr sf ; + <Jdd)l&tot conveys information regarding the 
onset of s-unitarity at very high energies. Assuming diffraction to be exclusively GW, we 
obtain[TT] Rd < 0.5. Multi IP induced diffracion is not included in Rd since it originates from 
G™ k . Hence their non screened high mass cross section is suppressed by its survival probability 
which decreases with energy. In GLMM Rd < 0.5, decreasing slowly. In KMR(07) Rd > 0.5, 
increasing slowly with energy. The partial information available on KMR(08) suggests that its 
Rd grows even faster. 
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6 Conclusions 



This presentation centered on the phenomenology implied by multi Pomeron dynamics incor- 
porated in soft diffraction and its consequences on soft scattering. 

The concept of IP enhancement was triggered, at the time, by the assessment that G^p is 
not too small. Our view of the Pomeron got more focused with the updated data analysis of 
soft scattering in which we get Ajp ~ 0.30 — 0.35 and a' P ~ 0.01. The implied KMR Pomeron 
is hard enough to be treated partonically, in which the traditional classification of the soft IP 
as a simple J-pole and the hard IP as a branch cut in the J-plane is maintained. GLMM went 
further on identifying the soft IP with the hard IP. This is, clearly, a fundamental theoretical 
issue which should be further investigated. 

As it stands this dynamics is compatible with the data, but at present we can not support it 
with a decisive verification. A prediction shared by GLMM and KMR is the expected significant 
reduction, compared with non screened predictions, of a tot and ael at the LHC. In my opinion 
this measurment is of a critical value. 

Decisive experimental signatures of -P-enh are expected, essentially, above the Tevatron. 
Consequently, we should be prudent when evaluating phcnomenological models which reproduce 
the Tevatron data well. This is, obviously, required of a successful model, but is definitely not 
sufficient. 
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